Period doubling in the 2D Antiferromagnet. New gauge-fields and their anomalies. 
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We construct an effective gauge field theory model describing a 2D quantum antiferromagnet in the 
flux phase ground state. Due to period doubling the number of gauge fields is four rather than just 
one. These additional gauge fields correspond to field configurations in which the flux is staggered, 
which should be taken into account to enforce single occupancy within the magnetic unit cell. We 
calculate a new type of quantum anomaly, present in such a model. In particular this leads to a 
generalized local flux hypothesis, which is tested numerically. 



Although there is considerable evidence to believe that 
the ground state of the two dimensional spin-^ Heisen- 
berg antiferromagnet on a square lattice with nearest 
neighbor interactions has Neel order, the ground state of 
the doped antiferromagnet is still very far from being un- 
derstood. We believe that one of the possible universality 
classes is correctly described by the so-called flux phase 
which is characterized by non-zero chirality of spins. Mo- 
tivated by this, ane can apply the mean field approach, 
first proposed ir£l, to the 2D Heisenberg model. 



This yields the following effective Hamiltonian: 
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Here is to be found self-consistently in such a way 
that it minimizes the energy of the system. To exclude 
the extra non-physical degrees of freedom one has to im- 
pose an additional constraint 
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There are several knowjusolutions for A^ which give local 
minima of the energyBcl. The real question, however, is 
which of thoin will be preferred by doping. One possible 



candidate isEl Ai 



AAi 



• n 



plaquctte 



-1. This 



state is called the flux phase because the effective Hamil- 
tonian (|^) — without constraint (||) — describes fermions 
in a background magnetic field, with flux ir (in units 
where the elementary flux quantum — is 27r) per pla- 
quette. We will focus on this state in the present paper. 
One can easily find the spectrum of such fermions to be: 
e(k) — ±2A-^/cos 2 k x + cos 2 k y . The filling factor — the 
ratio of the total number of electrons with given spin to 
the number of lattice sites — is 1/2. This means that in 
the mean field approximation the lower branch of this 
spectrum is completely filled with electrons of each spin. 
This mean field solution is not a very good approxima- 
tion to (|l|) because it does not take into account the local 
constraint (|^). This is why in this mean field the anti- 
ferromagnet looks like a metal while after taking into 



account the constraint it becomes an insulator. The way 
to deal with the constraint is to consider the fluctuations 
around the saddle-point mean field with 7r-fluxHLj. More 
precisely these are fluctuations of the phase of A, or Aij, 
that restore the single occupancy constraint (^). One can 
easily see this, if one notices that the fluctuating gauge 
field forces the current to be zero on every lattice bond. 
In what follows we restrict ourselves to phase fluctuations 
and adopt the common point of view that fluctuations of 
the modulus of Hubbard-Stratanovich field |A| are not 
relevant for the low energy sector of the theory because 
they have a gap with size of order J, see egjj. 

From now on we will leave spin out of the discussion, 
since in the models that we will discuss, both up- and 
down-spin electrons are decoupled. Let us start from 
mean field theory (0) neglecting the local constraint (||) 
and fluctuations of the gauge field. To lift the ambi- 
guity of the sign of the flux (— 7r = 7r mod 27r), we 
allow in addition to horizontal and vertical, also diag- 
onal hops of amplitude ^\A\ in all four diagonal direc- 
tions, effectively distributing the flux thmugh every pla- 
quette over its two constituent triangleax Then in the 
gauge A% m = , Av mn = f , A™ = and 

A™ = (-l)"f , where NE and NW refer to "north-east" 
and "north-west" diagonal bonds, corresponding to the 
7r-flux per plaquette we have the following spectrum: 



e(k) = ±2 1 AU/sin 2 k x + sin 2 k y 



'■k y (4) 



Because the periodicity of the lattice is doubled, due to 
the presence of half a flux quantum per plaquette, the 
Brillouin zone is half the size compared to the zero mag- 



netic field case: 
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< | . As one can see 



from the spectrum, at half filling, all contributions to the 
low energy dynamics of the theory come from the vicinity 
of the "Dirac" -points (k x , k y ) — (0, 0), (tt, 0). Now we in- 
troduce a slightly different representation which is more 
useful for studying low energy dynamics. We choose a 
2x2 unit cell on the square lattice, reducing the Brillouin 
zone further to \k x \, \ k y \ < 5. This choice of unit cell is 
convenient because in the reduced Brillouin zone there is 
only one Dirac point at k — (0,0), albeit doubly degen- 
erate. The low energy Hamiltonian can be constructed 
by simple Taylor expansion in momentum around th, 
point. This "lattice adopted basis" was proposed ir 
The result is the massive Dirac Hamiltonian: 
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(5) 



Here a M , (3 with /j, = 1, 2 are 4x4 matrices obeying the 
Clifford algebra. The wave function has 4 components, 
corresponding to the 4 sites in the unit cell (or 2 branches 
x 2 Dirac points). The model in its current form de- 
scribes two species of Dirac fermions, in 2+1 dimensions. 

Since the unit cell is now 2 x 2, we introduce the fol- 
lowing convention for decomposing an object O living on 
the L x L lattice, into four others living on a lattice of 
cells: 
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<D(r) =^A (l) C> w (r) (6) 

= C (1) (r) + {-) x+ yO {2) {Y) + {-) x O (3) {r) + (-) v O {4) {r) 

A = ±1 is the generalization of the sign-function: it dis- 
tributes the +, — signs within the unit cell according to 
one of the particular "staggered-ness" patterns, labeled 
by A: 
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and 



The main question we are going to address in this pa- 
per is how to incorporate fluctuations of the (lattice) 
gauge field into the continuum approximation (ph. The 
naive answer would be just to apply the minimal sub- 
stitution to (||). We argue that this is not enough to 
make the theory fully gauge invariant and thus to sat- 
isfy the no-double-occupancy constraint. Due to period 
doubling the replacement p — > p eA in (^) describes 
only gauge field fluctuations which are smooth inside the 
unit cell, i.e. of type A^ in (|J). These particular fluctua- 
tions alone would only partially restore the constraint (|^), 
just enforcing the total number of electrons inside each 
unit cell to be 4. For the constraint to hold locally on 
each site we have to include staggered gauge field fluctua- 
tions, of types A' 2 ', A^ 3 ) and A^ 4 ). These staggered fie 
give rise to additional fields in the continuum theoryQ' 
A theory that leaves those extra fields out, fails to de- 
scribe intra-cell charge-transfer correctly. Only the exact 
treatment of the single occupancy constraint can cure the 
mean field theory to restore insulating properties within 
the unit cell. 

One quick way to find which terms we should add to 
Hamiltonian (H) , is to consider the following: on the lat- 
tice we can perform a local gauge transformation, which 
might be different on each site within the unit cell. In 
the continuum this means, that we must be able to apply 
gauge transformations (which should leave TL invariant) 
that not only depend on r, but are also different for the 
4 components of the wavefunction ipi(r), i = 1 . . .4. I.e. 
i\) is multiplied not by a phase, but by a unitary diago- 
nal matrix, which is no longer proportional to the unit 
matrix. Using the previously introduced notation, we 
write such a general gauge field transformation compat- 
ible with period doubling (and smooth from cell to cell, 
where the word "cell" refers to a block of 2 x 2 plaquettes) 
as tp — > exp(i Ylt=i ^i^)^- ^ 1S the 4x4 matrix- 
representation of A, so #W = diag(Ag ' , Aq\\ A^q, A^\') 



e.g. 5^ = diag(l, — 1, 1, — 1). Then the simplest, mini- 
mal way to complete Hamiltonian (^|) to make it gauge 
invariant under all 4 of these transformations is to in- 
troduce 4 gauge fields E, F, G, and K in the following 
way: 



H = 




-*M„ 



f'N 



(3m e 



-iN 



(7) 



where = F^ + 5^E^ + S^G^ + S^K^ N = 
8^G 2 +8^K 1 , provided E, F, G and K gauge-transform 
as follows: F^ — > F M + (f>\ . . . K M — > + 4 , jj, = 1, 2. 
The particular form of N is dictated by the fact that we 
do not allow flux fluctuations through triangles. 

One can easily see that the usual minimal substitution 
gauge field can be identified as A\ = d\F\, A 2 = d 2 F 2 . 
The other fields correspond to staggered configurations 
of the magnetic field on the lattice as follows: 



d x d y (F 1 - F 2 ) 

{Ei - E 2 ) 
d x (G t - G 2 ) 
d y {Ki-K 2 ) 



+$M(r) 
-$(2)( r ) 

-i$( 3 )(r) ' 
+ i$ (4) (r) 



(8) 



and <f>W(r) is one of the symmetry components of the 
flux through plaquette r via (^). The important point 
here is that it is not possible to expand Hamiltonian (Q) 
to linear (or even any higher but finite) order in gauge- 
fields othep-than F and conserve full gauge invariance of 
the theoryt_J. 

We would now like to see how the electronic density re- 
sponds to fluctuations of these new fields. As we shall see 
staggered magnetic fields cause staggered density fluctu- 
ations, so we define the field induced staggered density 



as 



P W (r) + ^%),(i) 



E 

V>:-E<0 



The ipE form an orthonormal basis of eigenstates in the 
single particle Hilbert space. In the groundstate all neg- 
ative energy states are occupied, as indicated in the sum. 
To calculate these densities, we can use the following 
identity!-]: 



P 



W(k) 



-nn 27r. 



d 2 P 

(2tt) 2 



Tr(p\6 {l 



-H 



z 2 +H- 



■\p-k) (9) 



To make further progress, we have to consider the limit 
where fluctuations are small. Then, to first order in fields, 
we obtain after some algebra: 



p«(r) = 
P (2) W = 
P (3) (r) = 
pW(r) = 



sgn(m) (1 



i^ W (%) 

O Z7T 



4tt w 



(10) 

(11) 
(12) 

(13) 
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Here we kept only the zeroth order in — and the low- 
est order in momenta. The fact that (|l0|)-(|r^) survive 
when m — > oo indicates their quantum anomalous ori- 
gin. Therefore all staggered densities are expressed in 
terms of derivatives of staggered fluxes on the square lat- 
tice. Notice the non-trivial prefactors in the field-density 
relationships. Eq. ( |f0| ) is widely known as a manifesta- 
tion of the Chern-Simons term in the effective action of 
the ordinary gauge- field: ^ A)V x A. Relations ( |i"i"fjl3| ) 
are new, and represent the generalization of the Chern- 
Simons term to the other gauge- fields. 

Now we present some numerical support for these pre- 
dictions to see how they hold up on the original lattice, 
with Mean Field Hamiltonian (g). Our simulations were 
carried out on lattices of L x x L y sites {L x ,L y < 50), 
in the presence of periodic boundary conditions (other 
boundary conditions do not change our results quali- 
tatively). We let the electrons propagate in a uniform 
background magnetic field of flux Tr/plaquctte in addition 
to a small modulation. The absolute value of the hop- 
ping amplitude is constant and non-zero only for nearest 
and next-nearest neighbors. We apply a small, circularly 
symmetric perturbation to the uniform background flux 
$(n) = 7r, such that the additional flux through the pla- 
quette n is given by: <5$(n) = 7r/(n)A^)(n). Here /(n) 
is chosen to be gaussian: 7rAexp(— n-n/£ 2 ) + const. The 
constant is chosen such that the total flux through the 
entire lattice is unchanged. After diagonalizing the (sin- 
gle particle) Hamiltonian wc fill half of the spectrum with 
electrons and measure (n) . Notice that pM is defined 
only one quarter of the lattice sites. 

In our simulations we pick the ratio of the amplitudes 
of diagonal to NN hops to lie around 0.2, making the 
underlying field theory massive. The mass must be such 
that the gap is much smaller than the bandwidth m <C 
I, or y <c I. Since the result of the staggered density 
response (|Io|-|i3|) is obtained as an expansion in — we 

really have a window for m: £ _1 <C j <C I. In that 
case £ 3> 1, so that the continuum approach is valid. 
At the same time, £ has to be smaller than the lattice 
size L to avoid finite size effects: 1 <C £ <C L. Finally 
we should keep the amplitude of the perturbation small 
enough to remain within the realm of applicability of the 
linear approximation and anomaly calculation: 5$ <C 

™v*f «(£) 2 - 

Fig. |] shows an example where we have chosen a gauge 
fluctuation of the type S^ 4 \ Considering the a priori 
strong restrictions on the parameters to satisfy just men- 
tioned inequalities, the numerical experiments confirm 
the anomaly predictions rather well. In fig. ^ we in- 
vestigate quantitatively the response to this perturba- 
tion. Equation (|f^) tells us that the density response 
should be the derivative of a gaussian of the same width 
as the flux perturbation: d x e~( r '£' ~ xe~^ r /^> . The 
data points of fig. |^ — which show a slice through the lat- 
tice at fixed y = L/2 — fit this well with an amplitude 



that deviates < 10% from the predicted value of (|i~3]). 
In fig. H we repeat this exercise for the density response 
to the simplest gauge fluctuation: the field F (also cho- 
sen to be gaussian). According to ( |i"o| ) the predicted 
response should again be gaussian, and the data points 
fit this well, with an amplitude that deviates < 1% from 
the predicted value. The E-Held (type S^) obeys (plj), 
but for system sizes that we considered the amplitude is 
~ 30% too small. 

In conclusion, we showed that to maintain the single 
occupancy constraint at zero doping, extra gauge fields 
should be added to the continuum field theory description 
of (H). These fields represent staggered flux fluctuations 
on the lattice and are relevant for low-energy dynamics 
because of period doubling due to the 7r-flux mean field. 
The gauge invariance of the lattice gauge field theory 
demands that new gauge fields enter into the continuum 
description (^) as non-linear terms. The presence of these 
fields leads to anomalous terms for induced staggered den- 
sities. The situation, like the one described here, hasJpeen 
encountered in gravitational anomaly calculationsE3. In 
fact our new gauge fields E, G, K play the role of tetrades 
and spin-connections in gravity, except for the fact that 
due to fermion species doubling on the lattice, new gauge 
fields correspond to generally covariant transformations 
which in addition mix these two species of fermions in 
isospace. The new fields and new quantum anomalies 
found, must be important for any calculation of quanti- 
ties which involve intra-cell charge-transfer. In fact these 
fluctuations of staggered flux restore the insulating be- 
havior of the spin system within the unit cell. As the 
fluctuations determine the holon-doublon wavefunction, 
one should be able_to probe them directly by high-energy 
Raman scatteringE3. 
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FIG. 2. Density p' 4 ' plotted along a line of constant y for 
the case of flux perturbation of type 5' 4 ' taken from the data 
of fig. |l|. The solid line is the predicted density response from 
( p^[ ) on this slice. There is a finite size error in the amplitude 
of about 10%, illustrated by the dashed curve. 
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FIG. 3. Density p' 1 ' (r) plotted along a line of constant y for 
the case of a flux perturbation of type A — 0.05, £ = 14, 
50 x 50 sites. The slice of the 2D data is indicated by the 
fat line in the inset. The solid line is the predicted density 
response from (|l(]) on this slice. The fit is exact to within 
discretization errors. 




FIG. 1. Example of the application of a flux perturbation, (^4 = 0.05, £ = 14, 50 x 50 sites) of type S^K Panel a) shows 
the total applied flux, resulting from turning on field ^^(r) with a gaussian profile, as shown in the inset. Panel b) shows the 
total density response, and the inset shows its p^ 4 '(r) component. Panel c) illustrates what we would expect for p*- 4 ' based on 
formula (Jl3|) . 
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